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ABSTRACT 
The  reduction  of  circuit  equations  to  normal  form  for  numerical 
integration  is  considered  for  the  general  circuit  where  excess  reactive 
elements,  all  types  of  linearly  dependent  sources,  and  nonlinear  dissipa- 
tive  and  reactive  elements  are  present.    For  the  linear  circuit,  necessary 
and  sufficient  conditions  for  the  existence  of  numerical  solutions  are 
considered  and  stated.    For  the  nonlinear  circuit,  reduction  to  normal 
form  is  not  always  possible.     Numerical  solution  is  shown  to  be  simpli- 
fied if  certain  a  priori  conditions  are  satisfied  in  formulating  the  original 
circuit  equations.    A  new  systematic  reduction  procedure  is  presented 
for  obtaining  the  normal  form  equations.    This  procedure  is  also  extended 
to  a  new  procedure  for  obtaining  transfer  and  immittance  functions  of  the 
linear  circuit  from  a  proper  tree  formulation . 
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I.     INTRODUCTION 
The  formulation  of  circuit  equations  for  their  numerical  solution  is  of 

current  interest.    The  initial  contribution  was  made  by  Bashkow,  who 

12    3 
formulated  the  concept  of  the  proper  tree.    '    '        This  was  followed  by 

the  work  of  Bryant,  and  Kuh  and  Rohrer  who  refined  the  techniques.    The 
object  of  their  procedure,  which  was  based  primarily  upon  topology, 
is  to  obtain  the  circuit  equations  in  matrix  form  so  that  one  has  a  first- 
order  differential  equation.    The  approach  uses  the  state-variable  tech- 
nique, which  has  dominated  the  literature.    However,  it  should  be  noted 
that  other  formulations  of  circuit  equations,  such  as  the  TRAC  Program, 
do  not  use  state  variables. 

Once  the  circuit  equations  have  been  formulated  by  the  proper-tree 
concept,  they  have  the  following  form 

Ax     +      Bw      +     Cz      =     Dx      f     Ew      +     Fz      +     Gu 
where  x  is  the  vector  of  voltages  and  currents  associated  with  inde- 
pendent energy-storing  elements,  w  is  the  vector  of  currents  and  voltages 
associated  with  dissipative  elements,  z  is  the  vector  of  currents  and 
voltages  associated  with  excess  reactive  elements,  u  is  the  vector  of 
independent  sources,  and  the  capital  letters  are  coefficient  matrices. 

For  numerical  solution  it  is  necessary  to  alter  the  form  to  the  normal 
form    ,   X  =  Ax  +  Bu  for  the  linear  case  and  x  =  f(x,u,t)   for  the  non- 
linear case,  so  that  numerical    integration  can  be  applied  directly. 

It  has  been  shown  that  the  circuit  matrix  equation,  based  on  the 
formulation  of  the  proper  tree,  gives  the  relation  of  the  currents  and 
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voltages  of  each  circuit  element  to  the  currents  and  voltages  of  all  other 
circuit  elements.     Of  this  set  of  circuit  variables,  currents  and  voltages 
of  the  circuit  elements,  there  is  a  minimal  set  called  states  that  are 
required  to  describe  the  circuit.    If  those  variables  which  do  not  belong 
to  the  minimal  set  can  be  eliminated  the  resulting  matrix  equation  will  be 
of  normal  form,  or  an  equation  that  can  be  readily  converted  to  normal 
form . 

Chapter  II  is  a  study  of  the  circuit  that  contains  linear  circuit 
elements  and  linearly  dependent  sources.     The  formation  of  the  proper 
tree  is  discussed.    A  new  systematic  technique  for  matrix  reduction 
that  eliminates  the  unwanted  circuit  variables  is  developed.     The  circuit 
that  contains  no  dependent  sources  has  a  matrix  equation  that  always 
reduces  to  normal  form  and  reveals  that  the  number  of  states  is  equal 
to  the  number  of  independent  energy-storing  elements  in  the  circuit. 
The  introduction  of  dependent  sources  alters  the  coefficient  matrices 
of  the  circuit  equation.    These  alterations  make  the  reduction  process 
conditional  and  these  conditions  are  discussed.    In  addition,  it  is  found 
that  the  introduction  of  dependent  sources  may  alter  the  number  of  states 
required  to  describe  the  circuit.     Those  reactive  variables  that  are  not 
destined  to  become  states  are  designated  as  surplus.    It  is  shown  that 
the  surplus  reactive  variables  can  always  be  eliminated  from  the 
equations . 

In  Chapter  III  the  reduction  techniques  developed  in  Chapter  II  are 
utilized  in  establishing  a  new  procedure  for  deriving  driving-point 
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impedance/admittance  functions  and  transfer  functions  of  the  linear  cir- 
cuit.   This  procedure  enables  a  transfer  immittance  to  be  calculated 
directly  from  the  proper  tree  formulation.     Derivatives  are  replaced  by 
their  frequency-domain  equivalents. 

In  Chapter  IV  nonlinear  dissipative  elements  are  included  with  the 
linear  circuit  elements.    It  is  shown  that  if  normal  form  state  equations 
are  to  be  derived  the  currents  and  voltages  of  the  nonlinear  dissipative 
elements  must  be  restricted  in  their  relationship  to  the  voltages  and 
currents  of  other  circuit  elements.    This  relationship  may  be  controlled 
by  a  priori  circuit  modifications  that  increase  the  number  of  states. 

Chapter  V  extends  the  discussion  to  nonlinear  reactive  elements. 
Two  new  theorems  relating  to  the  conditions  for  the  elimination  of 
dissipative  variables  are  stated.    It  is  also  revealed  that  the  surplus 
reactive  variables  of  the  general  nonlinear  circuit  must  be  associated 
linear  elements  and  that  these  currents/voltages  must  not  be  related  to 
the  currents/voltages  of  nonlinear  elements  if  normal  form  equations  are 
desired. 


II.     REDUCTION  OF  LINEAR  CIRCUIT  EQUATIONS 
TO  NORMAL  FORM 

State  equations  are  often  used  in  circuit  analysis.     Normal  form 

state  equations,  x  =   f(x,u,t),    may  be  obtained  from  the  matrix  equation 

of  the  circuit's  currents  and  voltages. 

A.      FORMATION  OF  THE  CIRCUIT  MATRIX  EQUATION 

In  order  to  derive  circuit  equations  in  a  consistent  manner  a  proper 
tree  is  used.'    Any  tree  for  the  circuit  may  be  used  to  obtain  the  circuit 
equations.     However,  if  other  than  the  proper  tree  is  used,  additional 
steps  of  transposing  the  variables  to  the  proper  side  of  the  equation  are 
necessary  in  the  reduction  process.    The  proper  tree  is  drawn  as  follows 
All  voltage  sources  are  designated  as  tree  branches  and  all  current 
sources  are  designated  as  links.    This  precludes  a  loop  of  voltage 
sources,  which  necessarily  must  be  constrained  so  that  the  loop  can  be 
eliminated,  or  a  cut-set  of  current  sources,  which  must  be  constrained 
so  that  the  cut-set  can  be  eliminated,  as  illustrated  in  Figure  1.    All 
capacitive  elements,  less  the  minimum  required  to  prevent  the  formation 
of  a  loop,  are  designated  as  tree  branches  (C);  those  exluded  become 
links  (elastance,  S) .     All  inductive  elements,  less  the  minimum  required 
to  prevent  the  formation  of  a  cut-set,  are  designated  as  links  (L);  those 
excluded  become  tree  branches   (inverse  inductance,  T  ).     Resistive 
elements  which  do  not  complete  loops  are  designated  as  tree  branches 
(conductance,  G);  the  excess  become  links   (resistance,  R) .    Thus  the 
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order  of  precedence  for  drawing  tree  branches  is  (1)  voltage  sources, 
(2)  capacitive  elements,   (3)  dissipative  elements,  and  (4)  inductive 
elements . 

It  is  desirable  to  give  consideration  to  elements  whose  currents 
and  voltages  are  used  as  controlling  terms  for  dependent  sources.    When 
a  source  is  voltage  dependent,  the  element  across  which  the  controlling 
voltage  is  developed  should  be  a  tree  branch  if  this  constraint  does  not 
upset  the  precedence  for  forming  a  proper  tree.    When  a  source  is  current 
controlled,  the  element  through  which  the  controlling  current  flows  should 
be  a  link.     Should  these  arrangements  not  be  possible,  the  controlling 
element's  branch  current  may  be  written  in  terms  of  the  link  currents  of 
elements  forming  a  cut-set  with  the  controlling  element.     Similarly,  a 
link  element's  voltage  may  be  written  in  terms  of  the  branch  voltages  of 
elements  forming  a  loop  with  the  controlling  element.    Should  the  net- 
work be  too  complex  to  readily  perform  these  substitutions  and  arrange- 
ments, this  special  consideration  should  be  disregarded.     Most  circuits 
have  very  few  dependent  sources  in  comparison  to  other  circuit  elements. 
This  consideration  causes  the  dependency  constants  to  appear  in  the 
circuit  equations  with  the  constants  dictated  by  the  network  topology. 

Kirchhoff's  voltage  law  is  written  for  each  fundamental  loop  defined 
by  links  L,  R,  and  S.    Kirchhoff's  current  law  is  written  for  each  funda- 
mental cut-set  defined  by  tree  branches  C,  G,  and  r.    AH  voltage  sources 
appear  in  loop  equations  and  all  current  sources  appear  in  cut-set 
equations.    Figures  2  and  3  illustrate  the  conceivable  loops  and  cut-sets. 
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A  loop  defined  by  an  L  may  contain  dependent  voltage  sources,  indepen- 
dent voltage  sources  ,  C's  ,  G's  and  T  's  .    A  loop  defined  by  an  R  may 
contain  all  of  the  foregoing  with  the  exception  of  F's.    A  loop  defined  by 
an  S  does  not  contain  G's  and  F's.    A  cut-set  defined  by  a  C  may  contain 
L's,  R's,  S's,  and  current  sources.     Cut-sets  defined  by  G's  do  not 
contain  S's  and  cut-sets  defined  by  F's  do  not  contain  S's  and  R's.     By 
using  a  proper  tree  other  loops  and  cut-sets  are  prohibited. 

In  the  equations  which  follow,  F  (with  appropriate  subscripts) 
relates  link  and  tree  branch  voltages  and  link  and  tree  branch  voltages 
and  link  and  tree  branch  currents  according  to  the  topology  of  the  circuit. 

A  prime  designates  the  transpose  of  a  matrix,     e    ,  e    ,  and  e     represent 

L        K  b 

independent  voltage  sources,    j     ,  j     ,    and  j      represent  independent 

F      G  C 

current  sources,    r  and  /i  (with  appropriate  subscripts)  are  submatrices  of 
the  proportionality  constants  of  linearly  dependent  voltage  sources,     h 
and  g  (with  appropriate  subscripts)  are  the  proportionality  constants  of 
the  linearly  dependent  current  sources,    h  and  ^  are  dimensionless  where- 
as r  and  g  have  dimensions  of  ohms  and  mhos  respectively.     In  the 
matrix  equation  which  follows,  fundamental  loop  voltage  equations  and 
fundamental  cut-set  current  equations  are  written  alternately  as  indicated. 
Thus,  the  circuit  equations  become  (1). 
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(1) 


B.      THE  REDUCTION  PROCESS 

The  reduction  process  includes  the  elimination  of  the  variables 
associated  with  the  dissipative  elements,  the  transposition  of  the 
variables  associated  with  the  excess  reactive  elements,  and  the  elimina- 
tion of  the  coefficient  matrix  from  the  left  side  of  the  reduced  equations  . 
The  order  of  the  first  two  steps  is  unimportant. 
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Partitioning  (1)  as  indicated  and  lettering 

p.] 


[x]  = 


rwi= 


V 


c 


R 


V 


G 


V, 


iZ    = 


(2) 


enables  (1)  to  be  written  as  (3). 
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Q  Q 

^12  ^13 
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^22  ^23 

Q  Q 

^32  ^23 


X 

e 

X 

w 

+ 

e 
w 

i 

e 
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(3) 


In  (3),  X  designates  the  state  variables,  w  designates  the  "unwanted" 
dissipative  variables,  and  z  designates  the  excess  reactive  element 
variables.    Appendix  A  relates  the  submatrices  of  (3)  to  the  submatrices 
of  (1). 

1 .     Elimination  of  Dissipative  Variables 

The  following  rearrangement  of  (3)  shows  the  relation  of  the 
variables,  w,  to  the  other  variables  and  independent  sources. 


f(Pl2-'^12'" 


!^^22      ^22^ 
1(^32  -•532> 


w       = 


Q       Q 
^21  ^23 

Q       Q 
^31  ^33 


X 
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^    p  1 
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21     23 

P       P 
31     33 
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re  ■ 

i    X 


;    w 
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!  e 
L  z  j 


(4) 


To  have  a  determinate  solution  for  the  variables,  w,  the  rank  of 
the  matrix  multiplying  w  in  (4)  must  be  equal  to  the  number  of  dissipative 
elements  in  the  circuit.    This  matrix  has  dimensions  n  by  m .    Where  n  is 
the  total  number  of  circuit  elements  and  m  is  the  total  number  of 
dissipative  elements  in  the  circuit.     Because  of  the  dependent  sources, 
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whose  proportionality  constants  may  be  any  real  number^  there  exists  the 
possibility  that  the  matrix  multiplying  w  in  (4)  will  have  rank,  less  than 
m.    In  this  event  the  solution  for  the  variables,  w,  will  be  indeterminate. 

Let  the  rank  of  the  matrix  multiplying  w  in  (4)  be  m-i  because  of 
row  dependence.    Rearranging  the  rows  of  (4)  to  group  the  rows  such 
that  the  first  tier  is  dependent  upon  the  second  tier  leads  to  the  following. 


^ii^i2:'"(i) 


"ll^2 


x 


1^21  ^22j:^™(m-i)J       [N^i  N^^^. 


11^12:^^1^      K 


T       T 
21     22 


1  ) 

z  i 


(4a) 


M^_  is  a  square  matrix  of  rank  m-i  and  dimensions  m-i  by  m-i.    The 


enclosed  subscripts  denote  the  dimensions  of  the  w  vector  partitions. 

Since  the  first  tier  of  M  is  a  linear  combination  of  the  second  tier,  there 

must  exist  some  matrix  K  such  that  K     M^^  M^^      =   Tm,  ,  M,^     . 

L   22       22  J        L     1  i       12  _ 

Multiplying  the  second  tier  of  (4a)  by  -K  and  adding  the  results  to  the 


first  tier  yields 
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(4b) 


Th8  second  tier  of  (4b)  yields  the  indeterminate  solution  for  w. 
™(m-i)=  ^2"2('^21^^N224-T2iA-  T^^z   -  M^^w^.,   +  e^)  (4o) 

Recalling  that  there  were  n  elements  in  the  circuit,  m  of  which 
were  dissipative,  it  can  be  seen  that  the  first  tier  of  (4b)  contains  n-m+i 
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equations  relating  to  n-m  reactive  variables.     There  is  no  guarantee  that 
solutions  obtained  from  a  set  of  n-m  of  these  equations  would  be  con- 
sistent with  solutions  obtained  from  a  different  set.    Assuming  that  the 
state  equations  obtained  from  the  n-m+i  equations  were  consistent,  the 
substitution  of  these  states  and  their  derivatives  into  (4c)  does  not  yield 
a  determinate  solution  for  the  dissipative  variables.    Therefore,  it  is 
necessary  that  the  coefficient  matrix  of  the  w  variables  in  (4)  have  rank 
equal  to  the  number  of  dissipative  elements  in  the  circuit  if  the  w 
variables  are  to  be  eliminated. 


The  submatrix  |  P     -Q        |  has  dimensions  m  by  m  .     If  the  rank  is  m, 


then  the  solution  for  w  is  found  in  the  second  tier  of  (4). 

"  =  <P22-^22'''  '^21='  "  ^23*  -^21*"   ^23^   ^  «w' 
Substituting  (5)  into  the  first  and  third  tiers  of  (4)  yields 
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(7a) 
(7b) 
(7c) 
(7d) 
(7e) 
(7f) 
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^     =     P33  ^<Q32-P32'  <P22-^22'"'  ^23 
^1      =     ^x     "<^12-Pl2'  (^22-^22'"'  «w 
^2     =    \    ^"^32-^32'  <P22-^22'"'  ^w 


(7g) 
(7h) 

(71) 
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In  the  event  that  (P     -Q      )   is  singular,  the  matrix  inverse  used 
in  (5)  would  not  exist.    However,  the  n  rows  of  (4)  may  be  rearrange  to 
provide  a  square  matrix  of  rank  m.    This  matrix  may  then  be  used  in  lieu 
of  (P     -Q    „)   in  an  equation  similar  to  (5)  yielding  results  similar  to 
(6).    Failure  to  achieve  a  suitable  arrangement  of  the  rows  (4)  leads  to 
an  indeterminate  solution  for  the  variables,  w.    Since  these  variables 
are  resistor  currents  and  voltages  and  must  remain  finite,  indeterminate 
solutions  are  not  anticipated  in  circuit  analysis  problems. 
2  .    Transposition  of  Excess  Reactive  Variables 

To  achieve  the  normal  form  equations  from  (6),  z  and  z  must  be 
transposed.    This  transposition  yields 
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(8) 


J  L  ■  2. 

The  variables,  z,  are  currents  and  voltages  of  reactive  elements 
that  were  deemed  excess  in  the  formation  of  the  proper  tree.    In  a  circuit 
free  of  dependent  sources  these  variables  would  never  become  state 
variables.    Because  of  the  dependent  sources  (8)  shows  that  this 
distinction  can  no  longer  be  made.    Thus,  if 
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X 


A    = 


L^J 


-77       6 


;    B  = 


r^  ci 


u   = 


V, 


J 


L    ^J 


(9) 


(8)  becomes 

Ay    =    By    +    U  (10) 

Further  reduction  depends  upon  the  rank  of  the  matrix  A.     If  A 
is  nonsingular,  then  (10)  yields  the  normal  form 


-1  -1 

y    =    A      By    +    A      U 


(11) 


If  the  matrix  A  is  singular,   (11)  is  not  valid  and  additional  steps  in  the 
reduction  process  are  required. 

3 .     Elimination  of  Surplus  Reactive  Variables 

The  matrix  A  in  (10)  can  be  singular  because  of  linear  dependence 
between  the  rows  or  between  the  columns  or  both.    If  the  rows  of  (10)  are 
regrouped  forming  one  tier  which  consists  of  the  independent  rows  of  the 
matrix  A  and  a  second  tier  with  the  remaining  rows  and  these  tiers  are 
partitioned  to  form  square  submatrices  along  the  principal  diagonals  of 
A  and  B,  the  following  results. 


A  A,,n 

ii      id   i 

di     dd  • 
L.  _i 


y.  ! 
1 


d 


11        id 

i^di  ^dd 


^1 

u 

i 

^d 

+ 

"d 

(12) 


In  (12)  the  subscripts  i  and  d  denote  independence  and  depen- 
dence respectively. 

Since  the  second  tier  of  the  matrix  A  in  (12)  is  a  linear 
combination  of  the  first  tier,  there  exists  some  matrix  K  such  that 

.    Multiplying  the  first  tier  by  -K  and  adding  the 


[a      A,,^,=KrA..  A.^ 
L  di     ddj       L  11     id_ 


18 


results  to  the  second  tier  yields 

B 

id   :  i'i  i  ^ 

0      0       ^  ^y 


^A..   A.    -Ty.l 
;     11      id         1 


11 


id 


(B,.-KB,.)         (B,,-KB.  ,) 
di        11  dd        id 


y. 
1 

+ 

u, 

1 

^d^ 

_"d_ 

(13) 


Assuming  that  (B  ,   -KB.  J   is  nonsingular,  the  second  tier  of  (13)  yields 
dd        id 


y^      =      (B,,-KB.  ,)    ^      (KB..-B,,)y.      +     KU.      -      U, 
■^d  dd        id  11      di     1  1  d 


(14) 


Substituting  (14)  and  the  assoicated  derivative  equations  into  the  first 
tier  of  (13)  yields 


A*y       =      B*y       +     C*U       +     D*U,      +     E*U.      +     F*U^ 
i  i  i  did 


(15) 


where 

A* 

B* 
C* 


A..      +     A.  .(B,  ,-KB.  ,)  (KB..-B,,) 
11               id     dd        id  11      di 

B..      +      B,  ,(B,,-KB.  ,)"-^  (KB..-B,.) 
11               id     dd        id  11      di 

B,  ,  (B,,-KB.  ,)""^K      +     I 
id     dd        id 


D*  = 
E*  = 
p*      = 


-B.  ,(B,,-KB.  ,) 
id     dd        id 


-1 


-A.  ,(B,,-KB.  ,)"^K 
id     dd        id 


A,  ,(B,,-KB.  ,) 
id     dd        id 


-1 


(16a) 
(16b) 
(16c) 
(16d) 
(16e) 
(16f) 


If  A*  is  singular  because  of  row  dependence  the  process  used  to 
derive  (15)  from  (10)  is  repeated.    Each  repetition  of  this  process  may 
yield  higher-order  derivatives  of  the  source  terms.    If  A*  is  singular 
because  of  column  dependence  a  process  similar  to  that  used  to  derive 
(15)  is  employed. 

Assume  that  the  matrix  A  in  (10)  is  singular  because  of  depen- 
dence between  the  columns.    Regrouping  and  partitioning  (10) 
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accordingly  yield  an  equation  identical  to  (12)  where 

A. 


id 


A 


dd 


!     11 


A 


di 


H 


(17) 


Multiplying  the  first  tier  of  (12),  partitioned  according    to  column 


dependence,  by  A  . .     yields 


K 


^di     ^dd 


■J        


A   ^B  A   ^B    , 

ii     ii  ii     id 


1 


B 


di 


B 


dd 


y. 


n 


+ 


11  1 

u. 


(18) 


Multiplying  iihe  first  tier  of  (18)  by  -A       and  adding  the  results  to  the 


second  tier  yields  (19). 


_°  °J  '1^ 


11    11 


a".-^b.^ 

11     id 


,y. 


(B,.-A    A,/b..)  (B,,-A    A./b.  ,)|    ;y 
di      di   11     11        dd     di   11     id      '_  c 


+ 


A-\. 

11  1 

u^-a^.a"/u. 

d      di   11      1 


(19) 


Substituting  the  solution  for  y     from  the  second  tier  of  (19) 


and  the  associated  derivative  equations  into  the  first  tier  yields 
*Ay.      =      *By.      +      *CU.      +      *DU,      +      *EU.      +      *FU 


d       (20) 


where 


*A        =     I   -  K(B^^-A^.A./b.^)    ^    (B^.-A^.A./b..) 
dd      di   11     id  di      di   ii     ii 

*B        =     AT.'^B..-A~.-^B.^(B^^-A^.Ar/B.^)~'^   (B^.-A^.aTI^B.  .) 
11    11      11     id     dd      di   11     id  di      di    ii    ii 

*C       =     a"/    +    Ar.'^B    ,(B,,-A^.A".-^B    J""^  A^.A"-^ 
11  11     id     dd      di   11     id  di   ii 

*D       =      -a"/b.  ,(B,,-A^.a"/b.^)"^ 
11     id     dd      di   11     id 

*E        =      -K(B,,-A,.a""^B.  ,)"-^A^.Ar,-^ 
dd     di   11     id         di   11 

*F        =     K(B,,-A,.a".-^B.  ,)~-^ 
dd      di  11     id 


(21a) 
(21b) 
(21c) 
(2  Id) 
(21e) 
(2  If) 
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If  the  matrix  *A  is  singular  and  has  row  dependence,  the 
process  beginning  with  (12)  is  repeated.     If  there  is  only  column  depen- 
dence, the  latter  part  of  the  process  is  repeated. 

In  the  foregoing  it  was  assumed  that  the  matrices    (B    ,-KB     ), 

used  in  (13)  through  (16),  and  (B      -A    A~.  B     ),  used  in  (19)  through 

(21)^  were  nonsingular.    If  these  matrices  are  singular,  there  will  exist 
some  rearrangement  of  the  columns  of  (13)  and  (19)  that  will  provide 
nonsingular  matrices  that  can  be  used  in  place  of  (B    ,-KB.  ,)   and 

(B  ,   -A, .A..  B.  ,).    This  statement  is  proved  by  contradictions  that  show 
dd     di   11     id 

that  the  rank  of  the  second-tier  submatrix  of  y  coefficients  in  (13)  or 

(19)  is  equal  to  the  number  of  variables  to  be  eliminated.    First,  it  is 

shown  that  there  are  no  dependent  rows;  secondly,  it  is  shown  that  the 

number  of  independent  columns  is  equal  to  or  greater  than  the  number  of 

variables  to  be  eliminated. 

Let  the  following  represent  the  arrangement  of  (13)  or  (19)  in 

which  the  independent  rows  of  the  second  tier  of  the  y  coefficient 

matrix  have  been  separated  from  the  dependent  rows,  the  dependent  rows 

being  placed  in  a  third  tier. 


(22) 


Since    Fb'' ,    B"      i  =  K  [b*     B*      '[  ,    K  being  the  matrix  that  shows 
L  di     dd  J  L  di     dd  J 

the  linear  combination  of  the  second  tier  of  (22)  that  results  in  the  third 


A.. 

'A.. 

B 

B    , 

11 

id 

^  •    ' — 

ii 

id 

0 

0 

1 

= 

B* 

B*  , 

• 

di 

dd 

^d 

0 

0 

i-  _ 

B" 

B"  , 

di 

dd 

u. 

'y,  ' 

1 

X 

Jd. 

+ 

u* 

d 

2Ji 


tier,  multiplying  the  second  tier  by  -K  and  adding  the  results  to  the  third 


tier  yields 

A..      A.^ 

B,, 

B.^  ' 

i  u. 

11        id  i 

•       

11 

id  '  , 1 

i 

: 

y. 

I 

0         0 

1    _ 

B* 

B*     : 

+     !    U* 

. 

di 

dd     i 

;    d 

[o      0  _, 

0 

0  J  N. 

L^d 

-KU* 
d_ 

(22a) 


It  can  be  seen  that  the  third  tier  of  (22a)  contains  only  equations 
of  constraint  on  independent  sources.    Since  these  sources  cannot  be 
constrained  and  remain  independent,  it  must  be  assumed  that  the  second 
tiers  of  the  y  coefficient  matrices  in  (13)  and  (19)  contain  no  dependent 
rows . 

Assume  that  the  second  tiers  of  (13)  and  (19)  do  not  contain  a 
sufficient  number  of  independent  columns  in  the  y  coefficient  matrices 
to  perform  the  required  substitutions.    Let  the  following  represent  the 
rearrangement  of  the  second  tiers  of  (13)  or  (19)  which  provides  one 
square,  nonsingular  submatrix,  the  first  tiers  not  being  shown. 

(22  b) 

Let  C       be  the  nonsingular,  square  submatrix.     Since  the  second  set  of 
columns  of  C  is  a  linear  combination  of  the  first  set,  there  must  exist 
some  matrix  K  such  that 


1 

II 

^11 

^12" 

y* 

+ 

i  u* 

i 

0 

521 

^22.. 

y" 

L": 

c 


12 


22 


11 


C  L^ 

-1 


K 


(22  c) 


Multiplying  the  first  tier  of  (22b)  by  C  ^  ,     yields 


'22 


'"o^ 

'l 

K 

^y*" 

+ 

0 

_    _1 

_^21 

S2_ 

:  y"    1 

L  "■'. 

(2  2d) 

Multiplying  the  first  tier  of  (22d)  by  -C    ,  and  adding  the  results  to  the 

second  tier  yields 

K     '      Ty*    '  P      C"ju* 

(22e) 


-    -t 

0 

I 

1 

K 

W*' 

+ 

1 

0 

0 
L 

0 

y"    1 

L 

Again  constraint  equations  on  the  independent  sources  have  appeared. 
Thus,  it  must  be  assumed  that  there  are  a  sufficient  number  of  inde- 
pendent columns  in  the  second  tiers  of  the  y  coefficient  matrices  in  (13) 
and  (19)  to  perform  the  operations  required  to  eliminate  the  surplus 
reactive  variables.    Hence,  the  rank  of  the  second-tier  submatrix  of  y 
coefficients  in  (13)  and  (19)  is  equal  to  the  number  of  variables  to  be 
eliminated. 

C.     CONCLUSIONS 

In  a  circuit  without  dependent  sources  the  number  of  states  required 
to  describe  the  circuit  is  equal  to  the  number  of  capacitive  branches  plus 
the  number  of  inductive  links  in  the  proper  tree .    When  the  circuit 
contains  dependent  sources,  the  proper  tree  does  not  necessarily 
determine  the  number  of  states.    It  can  be  seen  by  comparing  (7)  to  (1) 
that  additional  states  maybe  required  if  the  dependent  sources  are 
controlled  by  currents  or  voltages  of  reactive  elements  made  excess  by 
the  formation  of  the  proper  tree. 
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Since  the  dependent  source  proportionality  constants  may  take  on 
any  value,  there  exists  the  possibility  that  fewer  states  than  predicted 
by  the  proper  tree  may  be  used  to  describe  the  circuit.     Because  of 
matrix  singularities  some  of  the  reactive  variables  become  surplus.    The 
removal  of  these  surplus  variables  introduces  derivatives  of  the  inde- 
pendent source  terms  and  reduces  the  number  of  states. 


D.     EXAMPLES 

1 .      Incr-ease  in  States  Due  to  Dependent  Sources 

Consider  the  circuit  of  Figure  4  which  illustrates  the  effect  of 
dependent  sources  that  requires  additional  states  to  describe  the  circuit. 
The  circuit  equations  for  Figure  4(a)  are  the  following. 


~1 


[h  hi 


L     t 

(2  ^L2 


^S^C3 


V 


SI 


I 


r3 


'J 


0       0       0        1  0        1 

0       0        11  0        1 

0-100  10 

-1-10       0  10 

0       0-1-1  0       0 

-1-10       0  0       0 


I 


LI 


L2 


V        i 
C2  • 


^C3! 


(23) 


^l^Sl 

Vr3 


With  the  absence  of  dissipative  elements  (22)  may  be  expressed  in  the 
form  of  (6)  in  which 


c  = 


LOO  0 

0       L^     0  0 

0       0       C^  0 

0       0        0  C 


;  a  = 


0 

0 

0 

1 

0 

0 

1 

i! 

0 

-1 

0 

Oi 

! 

-1 

-1 

0 

! 
0 

;  i8  = 


0     1 

0     1 

r^i  °] 

1     0 

_o    S. 

1  °J 

(23  a) 


M 


x= 


1 

0 

;    y  = 

0 

1_ 

;     C  =  ry   =  6    =    0 


0      0-1-1 

-1-1       0       0 

I _ 

Transposing  the  excess  reactive  variables,  the  equations 


(23  b) 


become 


L,     0       0        0  0        -L^ 

1  3 

0       L,     0        0  0        -L, 


Ll^ 


L2 


0       0       C^     0        -C^        0      !       V^^ 


0       0       0 


<=3     -^1       ° 


V 


0       0       0        0  0  0 

0       0       0        0  0  0 


C3 


SI 


0      0      0      10      0 

0      0      110      0 

0-10000 

-1-10      0      0      0 

0      0      1      1      1      0 


L. 


,       V„        '-1      0      0      0      0      1 


LI 


L2    I 


V 


C2 


V 


C3 


V 


SI 


I 


r3_ 

(24) 


(24)  is  in  the  form  of  (13)  and  the  variables  V       and  I        may 

be  eliminated.    With  the  dependent  sources  present  as  shown  in  Figure 
4(b). 


6    = 


rr     Ml 


h      g 


J 


C^        0 
0  L 


3J 


(25) 


and  the  equations  in  the  form  of  (24)  become 
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"] 


L,  0       0         0          0     -L, 
1 

:  0  L^    0         0          0     -L, 

I  2 

0  0      C^     0       -c,    0 


0     0     0       c^   -c^    0 

0       0       0         0       rC^     .uL^ 


L^ 


0      0      0        0      hC       gL^ 


LI 


L2 


V 


C2! 


V 


C3 


0  0  0  10  0 
0  0  110  0 
0      -10     0       0     0  I     V 


LI 


C2  1 


-1      -10     0      0     0  i  .  V 


C3 


V^  J      i   0        oil       1     0  1     V 


(2  6) 


i        I 


r3 


1         1     0     0       0     1  I     I 


SI 


r3 


If  6  is  nonsingular,  none  of  the  variables  can  be  eliminated.    If 


6   is  singular-  with  rank  of  one,  either  V    ,    or  I         may  be  eliminated. 

o  i  i  J 

Hence  it  can  be  seen  that  the  presence  of  the  dependent  sources  of 
Figure  4(b)  caused  an  increase  in  the  number  of  states  required  to 
describe  the  circuit. 

2.    Decrease  in  States  Due  to  Dependent  Sources 

Consider  the  circuit  of  Figure  5  which  illustrates  the  effect  of 
dependent  sources  that  reduces  the  number  of  states  by  causing  matrix 
singularities.    The  circuit     equations  for  Figure  5  are  the  following. 

r  Li. 
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L°2^GiJ 
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0 

-^ 

-1 

0 

0 

1 

0 

0 
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\ 

0    1 

J 

0 

0 

0 

0 

I, 


V 


c 


Rl 


V 


-1  L_ 


G2J 


[v] 


(27) 


The  equations  of  (27)  are  converted  to  the  form  of  (3).     Since 


there  are  no  excess  reactive  elements  P       =  Q^. 


^32  ^^32 


^33     =   ^^33     =     Pl3     =    '^13=P23  =  <=23     =    °- 
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[v]       (2  8) 


From  (28)  it  can  be  seen  that  (P     -Q      )   is  singular  and  can 

not  be  used  in  the  solution  for  the  dissipative  variables.     Noting  that 
P     -Q     )   is  nonsingular,  the  first  two  rows  of  (28)  are  used  to  solve  for 


these  variables . 


Rl 


V 


G2 


0  -1 

1  R 


1       '0       C        V^ 
L      L  -J       L  *^ 


I  0  c 

1-L        -CR 


V, 


(2  9) 


Substituting  (29)  into  the  second  tier  of  (28)  yields 


:G^       G^CRj 


L 


0        -1 
-1  0 


V, 


+ 


0 

L    J 


V 


(30) 


(30)  is  in  the  form  of  (13).    From  the  first  row  it  can  it  can  be 
seen  that  V      =  V.    Thus  substituting  V  for  V      in  the  second  row  of 

(30)  yields  the  final  result. 


°2"l     =      -\      -°2^°lV 


(31) 


^rom  (31)  it  can  be  seen  that  V      is  not  needed  to  describe  the 
circuit  of  Figure  5.    Had  the  proportionality  constant  r  been  any  other 
value,  V      would  have  been  a  state.    Thus  it  can  be  seen  that  dependent 
sources  may  reduce  the  number  of  states  required  to  describe  the 
circuit- 
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III.     LINEAR  CIRCUIT  EQUATIONS 
REDUCED  TO  IMMITTANCE  FORMS 

Linear  circuit  equations  derived  by  the  proper-tree  concept  can  be 
reduced  to  impedance,  admittance,  and  transfer  functions  if  the  circuit 
equations  are  augmented  by  output  equations  .    When  dealing  with 
transfer/immittance  functions  the  Laplace  operator  with  zero  initial  con- 
ditions must  be  used  in  lieu  of  derivative  notation.     If  the  circuit  equa- 
tions are  x  =  Ax  +  Bu  and  the  output  equations  are  y  =  Cx  +  Du,  then  x 
may  be  found  from  the  circuit  equations  and  substituted  into  the  output 
equations  to  yield  Y(s)  =    [C(sl-A)      B  +  D]  U(s).    This  is  an  extended 
determination  of  a  transfer  immittance  which  reduces  to  the  usual  one 
when  there  is  a  single  input  and  a  single  output.    This  procedure  requires 
one  to  obtain  the  state  equations  and  output  equations  and  from  them 
derive  equations  that  yield  the  desired  functions.    The  procedure  develop- 
ed here  enables  one  to  derive  the  desired  functions  directly  from  (1)  by 
reduction  techniques. 

A.      DERIVATION  OF  DRIVING  POINT  ADMITTANCE/IMPEDANCE 

If  a  network  is  driven  at  one  set  of  terminals  by  a  voltage  source 
and  a  proper  tree  is  drawn  of  the  network,  only  one  fundamental  cut-set 
will  be  described  by  the  source.    The  sum  of  the  currents  of  the  link 
elements  forming  the  cut-set  with  the  source  is  the  current  passing 
through  the  source.     Thus,  the  circuit  equations  may  be  augmented  with 
an  output  equation  which  expresses  the  source  current  in  terms  of  link 
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currents  .     Solving  for  the  source  current  in  terms  of  the  source  voltage 
results  in  an  equation  of  the  following  form  . 

I(s)    =  Y(s)V(s)  (32) 

I(s)  is  the  source  current,  Y(s)  is  by  definition  the  input  admittance 
function,  and  V(s)  is  the  source  voltage  function.     If  V(s)  is  unity  (32) 
becomes 

I(s)   =  Y(s)  (33) 

If  the  network  is  driven  by  a  unit  current  source,  only  one  funda- 
mental loop  will  be  defined  by  the  source  when  the  proper  tree  is  drawn. 
The  voltage  developed  by  the  source  is  the  sum  of  the  voltages  of  the 
branch  elements  forming  the  loop  with  the  source.     Thus,  the  circuit 
equations  may  be  augmented  with  an  output  equation  which  expresses  the 
source  voltage  in  terms  of  the  branch  voltages.    Solving  for  the  source 
voltage  in  terms  of  the  unit  current  source  results  in  an  equation  of  the 
following  form . 

V(s)      =     Z(s)  (34) 

V(s)  is  the  source  voltage  function,  and  Z(s)  is  the  input  impedance 
function. 

Writing  the  circuit  equations  in  the  form  of  (1),  with  LI       replaced 
by  sLI    ,  etc.  ,  and  performing  the  substitutions  given  by  (2)  yields  an 
equation  similar  to  (3),  as  follows. 
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Augmenting  (35)  with  the  output  equation  yields 
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where  6  is  the  output  current/voltage  to  be  determined  and  F      ,  F       , 

ox  OW 

and  F        are  terms  which  relate  the  output  to  the  circuit  voltages  and 

OZ 

currents.    Th^  coupling  terms,  F*  and  F**,  may  include  dependent 
source  constants  as  well  as  topological  factors. 

In  the  reduction  process  of  Chapter  II,  to  obtain  the  normal  form 
equations,  all  dissipative  currents  and  voltages  were  eliminated.    The 
elimination  process  caused  the  dissipative  parameters  to  appear  in  the 
coefficient  matrices.    The  output  equation  identifies  the  currents  and 
voltages  of  interest  and  to  obtain  an  immittance  the  proper-tree  equations 
should  be  reduced  by  eliminating  those  currents  and  voltages  that  do  not 
appear  in  the  output  equation.     Thus,  after  writing  the  equations  in  the 
form  of  (36),  the  unwanted  variables  are  eliminated.    The  reduction 
process  for  these  variables  is  identical  to  the  process  used  to  eliminate 
the  dissipative  variables  in  the  derivation  of  normal  form  equations, 
(see  Appendix  C) 

The  elimination  of  the  unwanted  variables,  those  for  which  the  total 
output  coupling  (F**  -F*)  is  zero,  reduces  the  dimensions  of  the  matrix 
equation.    In  addition,  the  reduction  may  cause  the  coefficients  of  the 
reduced  equation  to  become  polynomials  or  ratios  of  polynomials  of  the 
variables.    The  reduced  equation  is  of  the  following  form. 
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where  the  subscript  r  denotes  the  alteration  of  the  coefficient  matrices 

and  y  denotes  the  circuit  variables  that  were  not  eliminated.    Any  of  the 

elements  of  P  .  Q   ,  F*.  and  F**  may  now  be  functions  of  the  variable  s. 
r       r       r  r 

However,  this  notation  has  been  omitted  for  simplicity  in  writing. 
Solving^for  y  in  the  first  tier  of  (36a)  yields 
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M    =    [V^r]         k] 


(36b) 


Substituting  this  solution  for  y  into  the  second  tier  of  (36a)  and  solving 
for  6  yields 

e=  rp**-p*nrp-Q  i"^e  (37) 

Lr         rJLr      rj        y 
If  9  is  the  current  of  a  unit  voltage  source,  then  the  operations 
indicated  by  (3  7)  yield  an  admittance  function.    If  6  is  the  voltage 
developed  by  a  unit  current  source,  (37)  yields  an  impedance  function. 

B.      DERIVATION  OF  TRANSFER  FUNCTIONS 

Network  transfer  functions  may  also  be  derived  from  the  foregoing 
equations  by  letting       Q  become  a  voltage  or  current  at  terminals  not 
directly  connected  to  the  source.    Care  must  be  exercized  when  the  out- 
put is  to  be  an  open-circuit  voltage  as  in  the  case  of  a  transfer  impedance 
function.    There  may  be  elements  of  the  circuit  that  are  in  series  with  the 
open-circuit  terminals  that  do  not  affect  the  open-circuit  voltage.    These 


51/ 


elements  should  be  excluded  from  the  proper  tree  and  the  resulting 
circuit  equations. 

If  the  network  is  excited  by  a  unit  current  source  and  9  is  the  open-    , 
circuit  voltage  at  another  set  of  terminals,  then  (37)  yields  a  transfer 
impedance.    If  9  is  a  current,  the  right  side  of  (3  7)  becomes  a  current 
gain  function.     If  the  network  is  excited  by  a  unit  voltage  source  and  9 
is  a  current  through  another  set  of  terminals,  the  right  side  of  (37) 
becomes  a  transfer  admittance  function.    If  9  is  a  voltage,  a  voltage 
gain  function  is  derived. 

By  letting  0  become  an  output  vector,  both  input  and  transfer  func- 
tions may  be  derived  simultaneously.    Thus   (37)  may  be  extended  to 
multiple-port  circuits  with  multiple  excitations  by  identifying  the  sources. 
Since  the  output  equations  determine  the  unwanted  variables,  it  may  be 
necessary  to  perform  the  elimination  process  in  stages  to  obtain  various 
input  and  transfer  functions  . 

C .     SUMMARY 

The  steps  to  the  derivation  of  driving-point  and  transfer  functions 
are  summarized  as  follows: 

1.  determine  the  sources  needed  to  derive  the  desired  functions 

2.  eliminate  the  variables  that  appear  in  none  of  the  output 
equations 

3.  perform  the  operations  indicated  by  (37)  on  those  output  equa- 
tions where  the  output  is  coupled  to  all  the  remaining  variables 
including  those  variables  that  are  not  coupled  but  not  effective. 
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D.     EXAMPLES 


Consider  the  network  shown  in  Figure  6a.    Assume  that  the  input 
admittance,  Y,    ,  for  port  1-1'  is  desired.    Excite  the  circuit  with  a  unit 
voltage  source  at  port  1-1'.    The  proper  tree  is  shown  in  Figure  6b.    The 
augmented  circuit  equations  are 
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The  last  equation  is  the  output  equation  which  shows  the  funda- 
mental cut-set  relation  of  the  source  current  and  I  .  The  remaining 
circuit  variables  are  to  be  eliminated.  Elimination  of  the  dissipative 
variables  from  (38)  yields 
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Elimination  of  the  capacitive  variables  from  (39)  yields 
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Eliminating  I        and  performing  the  operations  indicated  by  (37) 


yields 
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where 
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Assume  that  the  two-port  3  parameter  matrix  for  the  network  of 
Figure  6a  is  desireci.  The  proper  tree  with  both  ports  being  excited  by 
current  sources  is  shown  in  Figure  6c.    The  augmented  circuit  equations 


are 
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From  the  output  equations  it  can  be  seen  that  I  _,  I   ,  ,  and  I   ^ 

LiZ        Rl  Rvj 

should  be  eliminated.    The  result  of  their  elimination  and  the  relocation 
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Performing  the  steps  necessary  to  obtain  the  form  of  (37)  yields 
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Completion  of  the  multiplication  indicated  by  (44)  reveals  the  desired 


parameters.  ,The  results  are 
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(46) 


The  solutions  have  been  left  in  a  form  that  allows  one  to  check  their 
accuracy  by  inspection. 

The  hybrid  parameters  could  be  derived  by  exciting  the  circuit  at 
one  port  with  a  voltage  source  and  a  current  source  at  the  other  port. 
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IV.     REDUCTION  TO  NORMAL  FORM 
WITH  NONLINEAR  DISSIPATORS 

Because  of  the  nonlinear  current-to-voltage  relations  of  semi- 
conductors, the  numerical  solutions  for  circuit  containing  these  devices 
may  become  complex  if  the  state  variable  approach  is  used.    The  circuit 
can  be  represented  by  a  set  of  normal-form  equations  that  are  augmented 
by  equations  which  also  must  be  solved  to  obtain  the  nonlinear  currents 
and  voltages  that  appear  in  the  normal-form  equations.    The  general 
form  of  these  equations  are 

X     =     Ax      +     B   u      +     Of    (v)      +      Df    (i) 

v      =     A^x      +     B^u      +     f3{i,v)  (47) 

i      =     A^x      +     B^u      +      f^a.v) 

where  the  f'x  are  nonlinear  functions  and  v  and  i  are  the  voltages  and 
currents  of  the  semiconductors.    If  f     and  f    are  zero,  then  (47)  reduces 
to  one  nonlinear  differential  matrix  equation  which  may  be  integrated 

directly. 

4 
Diodes  either  by  themselves,  or  as  part  of  the  Ebers-Moll 

equivalent  for  the  transistor  may  be  considered  in  terms  of  their  current/ 

voltage  characteristics.    For  the  solid-state  diode  it  is  possible  to 

write 
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The  exact  nonlinear  relation  is  unimportant  for  the  present  consideration 
so  that  the  form 

i,     =     f(v^)  (48a) 

d  d 

and  the  inverse  form 

V,     =     g(i^)  (48b) 

d  d 

are  used  here. 

Werther  and  Parker    and  Calahan     consider  the  diode  to  be  a 
current-dependent  voltage  source  or  a  voltage-dependent  current  source 
where  the  controlling  voltage  or  current  is  written  in  terms  of  the  voltages 
and  currents  of  other  circuit  elements  as  in  (47).    This  treatment  requires 
the  solution  of  network  equations  for  the  determination  of  the  controlling 
voltages  and  currents.    If  the  both  the  current  function,     (48a),  and  the 
voltage  function,   (48b),  are  single  valued  (the  tunnel  diode  is  exluded) , 
then  the  nonlinear  dissipative  element  may  be  treated  as  the  linear 
dissipative  elements  are  treated  when  writing  circuit  equations.    That 
is,  the  device  may  be  used  as  either  a  branch  element  or  a  link  element 
when  forming  a  proper  tree. 

A.      TOPOLOGICAL  CONSIDERATIONS 

The  nonlinear  elements  are  given  precedence  between  capacitive 
elements  and  linear  dissipative  elements  when  establishing  the  branches 
of  the  proper  tree.    When  writing  a  voltage  equation  for  a  loop  defined  by 
a  nonlinear,  dissipative  link  element,  a  term  on  the  left  side  of  the  equa- 
tion such  as  RT     is  replaced  by  the  voltage  v_,.    The  term  corresponding  to 
R  d 


I     of  the  right  side  of  the  equation  is  f (v  ) .     Similarly,  when  writing  a 
R  ci 

eufrent  equation  of  a  cut-set  defined  by  a  nonlinear,  dissipative  element, 

the  term  on  the  left  side  of  the  equation  is  i  ,  and  the  term  on  the  right  side 

d 

ii  g(i  ,).    The  matrix  equation  for  the  circuit  will  be  similar  to  that  of  (1) . 
d 


B.      ELIMINATION  OF  DISSIPATIVE  VARIABLES 

If  state  equations  are  to  be  obtained  the  dissipative  variables  must 
be  eliminated.    While  the  linear  dissipative  variables  may  readily  be 
eliminated,  the  elimination  of  the  nonlinear  variables  depends  upon  the 
relationship  of  the  currents/voltages  of  the  nonlinear  elements  to  the 
voltages  and  currents  of  other  circuit  elements  .    This  may  be  seen  by 
the  following  discussion. 

1 .    Nonlinear  Variables  As  General  Functions  of  Other  Variables 
Assume  that  the  nonlinear  voltages/currents  are  functions  of  all 
other  voltages  and  currents.    After  writing  the  circuit  equations  and 
completing  the  substitutions  given  by  (2)  the  equations  become 
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where  x,z,  and  w  are  given  by  (2)  and  P      and  Q.  ,  (i,j ,  =  1 ,2  ,3)  are 
given  by  Appendix  A.     P, ,   and  Q      (i  =  4,5;  j  =  1 ,2  ,3  ,4,  5  and  J  =4,5; 
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i  =  1^2,3,4^5)  are  matrices  which  relate  the  voltages  of  nonlinear 
dissipative  links  to  other  voltages  in  the  fundamental  loops  which  are 
defined  by  nonlinear  elements,  and  the  currents  of  nonlinear  branches  to 
the  currents  of  other  elements  which  are  in  the  cut-sets  that  are  defined 
by  nonlinear,  dissipative  branch  elements. 

The  first  two  tiers  and  corresponding  columns  of  (49)  may  be 
combined  by  transposing  z  and  z  .    Thus,  by  forming  y  from  x  and  z  as 
was  done  in  (9),   (49)  becomes  the  following. 
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As  was  discussed  in  Chapter  II,  the  w  variables  may  be  eliminat- 
ed if  the  rank  of  the  combined  coefficient  matrices  for  these  variables  is 
equal  to  the  number  of  these  variables.    Their  elimination  reduces  (50) 
to  the  following . 
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where,  if  (A     -B     )   is  nonsingular. 
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The  multiplication  of  (51)  by  [K]  ~      yields   (47).     Further 
reduction  depends  upon  the  elements  of  M      ,  M      ,  ^-39'  ^"^^  ^'^'^*    '^^^^ 
may  be  seen  by  writing  the  second  two  tiers  of  (51)  as  follows  . 

V,  K__    K_     ~^  K_  M  V  M        M  f(v   ) 

d      ^      22      23  _      21    ^-^    _^        21    ^-j  ^      2    ^        22       23  d     ^53^ 

^d  S2    S3  Si  ^31  ^3  ^32    ^33      ^^V 

Equation  (53)  suggest  several  undesirable  solutions  for  the 
nonlinear  currents  and  voltages.     Transcendental  solutions  may  be 
required.     If  4:hese  are  not  required,  solutions  which  reveal  that  a  non- 
linear current/voltage  is  a  nonlinear  function  of  other  currents/voltages 
may  result.    The  substitution  of  the  latter  solutions  into  the  first  tier 
of  (51)  yields  functions  of  functions.    If  M      =M      =M      =M       =  0, 

these  undesirable  solutions  cannot  exist.    By  comparing  (52)  with  (50), 
it  can  be  seen  that  the  desired  conditions  prevail  if  B„„=B     =B     =B     =B 

=B      =B     =B       =  0.    This  implies  that  the  fundamental  loops  defined 
by  nonlinear  elements  must  not  contain  other  dissipative  elements, 
either  linear  or  nonlinear.    Similarly,  cut-sets  defined  by  nonlinear 
elements  must  not  contain  other  dissipative  elements.    These  restric- 
tions must  also  apply  to  dependent  sources  which  are  controlled  by  the 
currents/voltages  of  dissipative  elements.    Thus,  if  desirable  equations 
are  to  result,  the  nonlinear,  dissipative  variables  must  be  restricted  to 
being  functions  of  reactive  variables  only. 
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2  .    Nonlinear  Variables  As  Functions  of  Reactive  Variables 

If  the  currents  and  voltages  of  nonlinear  elements  are  not  allowed 
to  be  controlling  terms  for  dependent  sources  and  these  currents/voltages 
can  only  be  related  to  capacitor  voltages  or  inductor  currents,  then  (50) 
becomes  the  following. 
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The  eliminations  of  the  w  variables  yields  the  following. 
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1^3. 


The  elimination  of  v  ,   and  i  ,   from  (55)   yields 

d  d 

K^^y  =  M^^y  +  V^    +  ^^^^(M^^y  +  V^)  +  M^3g(M3^y  +  V3) 


(56) 


If  K       in  (56)  is  nonsingular,  then  (56)  readily  yields  the  normal 
form  equations.    If  K, ,    is  singular  some  of  the  variables  in  y  must  be 
eliminated.    The  elimination  of  the  surplus  reactive  variables  may  be 
possible  but  very  difficult.    Assume  that  the  left  side  of  (5  6)  is  zero 
and  that  one  of  the  variables  y  is  to  be  written  in  terms  of  the  other 
variables  and  sources.    It  can  be  seen  that  the  solution  of  a  transcen- 
dental equation  may  be  involved.    To  investigate  the  conditions  neces- 
sary for  reduction  when  K       is  singular;  a  circuit  that  is  free  of 
dependent  sources  is  considered. 
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For  a  circuit  that  does  not  contain  dependent  sources  the  circuit 
equations  reduce  to  the  following. 
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The  matrices  Q      ,  Q      ,  Q      ,    and  Q        will  be  other  than  zero 
'  34         35         43  53 

whenever  a  fundamental  cut-set  is  formed  linear  dissipative  elements 

and  nonlinear  dissipative  elements.    Likewise,  Q       and  Q    .   will  be 

other  than  zero  whenever  a  fundamental  cut-set  is  formed  with  nonlinear 

dissipative  elements.    The  elimination  of  the  linear  dissipative  variables 

gives  rise  to  coefficients  for  f(v  )  and  g(i  ,)   in  the  reduced  equations 

causing  the  solutions  for  v  ,  and  i  ,  to  become  transcendental.    If  i  ,   is 

d  d  d 

eliminated  from  (57),  then  Q        and  Q        cause  the  solution  for  v     to 
become  transcendental.    Thus,  normal-form  state  equations  will  not 
result  from  the  reduction  of  circuit  equations  when  fundamental  cut-sets 
are  formed  with  nonlinear  dissipative  elements  and  other  dissipative 
elements . 

If  the  matrices  mentioned  in  the  preceeding  paragraph  were  zero, 
the  nonlinear  voltages/currents  would  be  functions  of  the  state  variables 
and  sources  only  and  all  dissipative  variables  could  be  eliminated.    The 
resulting  normal  form  equations  would  be  as  follows. 

X  =  Ax  +  BU  +  CU  +   Df(Ex  +  FU)   +  Gg(Hx  +  JU)  (5  8) 
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C.  CONCLUSIONS 

The  circuit  equations  for  circuits  containing  nonlinear  dissipative 
elements  do  not  readily  reduce  to  normal-form  state  equations  .     By  a 
priori  circuit  modifications  of  placing  an  inductor  in  series  with  those 
nonlinear  elements  which  have  loop  equations  formed  with  variables  that 
are  not  destined  to  become  state  variables  and  placing  a  capacitor  in 
parallel  with  those  nonlinear  elements  which  have  cut-set  equations 
formed  with  variables  that  are  not  destined  to  become  state  variables, 
the  nonlinear  elements'  voltages  and  currents  are  forced  to  become  state 
variables.    While  these  modifications  do  give  rise  to  additional  states 
they  allow  the  circuit  equations  to  be  reduced  to  normal-form  equations  . 

D.  EXAMPLE 

Consider  the  circuit  shown  in  Figure  7a.     The  circuit  equations  are 
the  following  . 
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This  equation  cannot  be  put  in  normal  form  because  of  the  coupling 
between  diodes.  If  the  circuit  is  modified  as  in  Figure  7b,  the  circuit 
equations  become 
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The  diodes  are  now  uncoupled  and  (60)  may  be  placed  in  normal 


form 
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(61)  is  obtained  by  eliminating  v      from  (60).     Normal-form  equations 


are  obtained  by  multiplying  (61)  by 

coo" 


0         C^        0 


0         0 


-1 


(62) 
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V.     REDUCTION  TO  NORMAL  FORM 
WITH  GENERAL  NONLINEAR  ELEMENTS 


When  nonlinear  inductors  and  capacitors  are  present,  the  circuit 
equations  may  be  reduced  to  the  following  form . 
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where 


X    = 


^, 
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;    £jx)      = 


V    \ 


-I 


{64a, b) 


and 


z  = 


Q, 


^, 


;    f,(z)    = 


V, 


^s'^s' 
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(65a, b) 


13',    is  the  flux  associated  with  inductors  and  Q     is  the  charge 
L  c 

associated  with  capacitors.    0     and  Q     are  functionally  related  to  the 

J-i  o 

inductor  currents  and  capacitor  voltages  as  given  by  (64b).    For  the 
linear  case  0     =  LI     and  Q   =CV   .     (65a, b)  give  the  general  relationships 
for  excess  reactive  elements. 

If  X  and  z  in  (49)  are  replaced  by  f   (x)   and  f   (z)  respectively,  the 
result  is  (63).    Since  z  and  f   (z)  may  be  transposed,  it  is  possible  to 
form 


^7 


y    = 


X 


and     f   (y) 

y 


(66a, b) 


If  f   (y)   is  used  in  place  of  y  in  the  equations  that  begin  with  (50),  the 
arguments  for  reduction  with  nonlinear  inductors  and  capacitors  would 
be  obtained.    Thus,  if  the  dissipative  variables  are  to  be  eliminated, 
the  conditions  for  the  elimination  of  these  variables  from  (49)  must  also 
apply  to  (63) . 

The  conditions  for  the  elimination  of  the  dissipative  variables  are 
stated  in  the  following  theorems. 

Theorem  1 .    If  a  coefficient  submatrix  for  the  linear  dissipative  variables 
in  the  equations  which  relate  these  variables  to  the  other  circuit 
variables  and  sources  has  rank  equal  to  the  number  of  linear  dissipative 
elements,  the  linear  dissipative  variables  may  be  eliminated  and  normal- 
form  equations  may  be  derived.     (For  proof  see  Chapter  IIBl) 

Theorem  2 .    If  each  of  the  voltages  or  currents  of  nonlinear  dissipative 
elements,  nonlinear  dissipative  variables,  are  linearly  related  to  only 
independent  sources  and  the  reactive  voltages  and  currents  that  are 
functions  of  variables  that  are  destined  to  become  state  variables,  then 
the  nonlinear  dissipative  variables  may  be  eliminated  and  normal  form 
equations  may  be  derived.     (For  proof  see  Chapter  IVB2) 

Having  eliminated  the  dissipative  variables  from  (63),  the  resulting 
equations  are  as  follows  . 

Ay  =   Bf    (y)      +      U      +     Cf^(Df   (y)      +     U)   +  Eg^(Hf   (y)    +  U)      (67) 
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If  A  in  (67)  is  singular  there  exist  surplus  states  to  be  eliminated. 

Since  these  states  do  not  appear  in  functions  f  ,   and  g  ,   which  result 

d  d 

from  the  elimination  of  the  nonlinear  dissipative  variables,  an  equation 
of  the  following  form  must  be  solved  to  eliminate  the  surplus  state 
functions  f   (y)   and  the  surplus  state  derivatives. 


-\2 


0       =     M      N 


y         Si 


where  the  subscripts  i  and  s  denote  the  state  functions  to  be  retained 

and  the  surplus  state  functions  respectively.    Wliile  the  surplus  state 

functions  f   (y)      may  be  written  in  terms  of  the  other  state  functions  and 
y      s 

sources,  the  solution  for  the  surplus  state  derivatives  involves  partial 
derivatives  of  both  the  surplus  functions  and  the  functions  to  be  retained, 

^^y^^^s  -1  ''^^^\  -1- 

-     -^ — ^     y       =     N   'm    '~^—  y.    +  N   ^U  (69) 

cy  s  Syi 

s  "^1 

From  (69)  it  can  be  seen  that  if  the  surplus  states  are  to  be  elimin- 
ated the  surplus  state  functions  must  be  linear  and  must  be  related  to 
retained  functions  that  are  also  linear.    Otherwise,  the  partial  deriva- 
tives must  be  entered  into  the  coefficient  matrix  of  the  retained 
derivatives  and  the  resulting  equations  become  nonlinear  differential 
equations  that  cannot  be  placed  in  normal  form. 
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APPENDIX  A 


General  Submatrix  Relations  to  Circuit  Equation  Submatrices 
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APPENDIX  B 
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(a)   Voltage  Source  Loops 
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(b)  Current  Source  Cut— Sets 

FIGI.       ELIMINATING  VOLTAGE  SOURCE  LOOPS 
AND  CURRENT  SOURCE  CUT-SETS 
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FIG  3.  CONCEIVABLE  CUT-SETS 
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(a)  Original  Circuit 
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FIG  4.       EXAMPLE  I.    INTRODUCTION  OF  STATES 
BY  MUTUALLY  COUPLED  EXCESS 
REACTIVE  ELEMENTS 
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FIG  5.       EXAMPLE  2.    REDUCTION  OF  STATES 
DUE  TO  DEPENDENT  SOURCE 
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FIG  6.   EXAMPLES  3  a  4.  IMMITTANCE  FORMS 
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FIG  7.  EXAMPLE  5.    INTRODUCTION  OF  ADDITIONAL 
STATE  TO  OBTAIN  NORMAL  FORM  REDUCTION 
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APPENDIX  C 

Elimination  of  Unwanted  Variables 

Consider  the  following  equation  and  the  elimination  of  the  variable 


w. 
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If  (A     -B„„)    is  nonsingular,  the  second  equation  may  be  written  as 

-  =    '^22-^22'"'   (^21^  ^  '*23^^<^)  -'^21^xM    "^3^  "  ''w' 

Substituting  this  value  for  w  into  the  first  and  third  equations  and 
combining  terms  yields 


^r^2-'^12"^22-B22'''  '^21;    '^13-^2-Bl2'<'^22-B22)'^'^23 

^r^2-'^32><'^22-^22'"'  ^21'    ^3-<^2-''32' <*22-B22'"' '^23_ 

Bu-^2-''l2)''^22-''22'"'  ''21;  ''l3-<^2-Bl2'<'^22-'^22'"'  ^2! 

^3r<^2-''32'"^22-B22'"'  '^21^  ^33"  ^2-832' ^2-^22'"'  ^23 


-1 


F^(x)l 


X 


^^(z) 


+ 


e   -(A     -B     )(A     -B_J       e 
X        12       12       22      22  w 

e    -(A      -B      )(A      -B.J"-^  e 
z        32      32       ZZ      22  w 


The  pattern  for  the  multiplication  and  combining  of  terms  is  readily 
established. 
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